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ABSTRACT 

We present simultaneous multicolor infrared and optical photometry of the black hole X-ray transient 
XTE Jl 118+480 during its short 2005 January outburst, supported by simultaneous X-ray observa- 
tions. The variability is dominated by short timescales, ~ 10 s, although a weak superhump also 
appears to be present in the optical. The optical rapid variations, at least, are well correlated with 
those in X-rays. Infrared JHK S photometry, as in the previous outburst, exhibits especially large am- 
plitude variability. The spectral energy distribution (SED) of the variable infrared component can be 
fitted with a power-law of slope a = —0.78 ± 0.07 where f„ tx There is no compelling evidence for 
evolution in the slope over five nights, during which time the source brightness decayed along almost 
the same track as seen in variations within the nights. We conclude that both short-term variability, 
and longer timescale fading, are dominated by a single component of constant spectral shape. We 
cannot fit the SED of the IR variability with a credible thermal component, either optically thick or 
thin. This IR SED is, however, approximately consistent with optically thin synchrotron emission 
from a jet. These observations therefore provide indirect evidence to support jet-dominated models 
for XTE Jl 118+480 and also provide a direct measurement of the slope of the optically thin emission 
which is impossible based on the average spectral energy distribution alone. 
Subject headings: accretion, accretion disks — binaries: close — stars: individual: XTE J1118+480 



1. INTRODUCTION 

The black hole X-ray transient (BHXRT) 
XTEJ1118+480 has proven to be a crucial object 
for understanding the class. It is one of the best objects 
for multiwavelength study as it suffers an exception- 
ally low interstellar extinction, allowing broad-band 
observations encompassing even the extreme-UV region 
lIHvnes et a,l.ll2000r iMcClrntock et a.ljl200ll IChatv et a.U 
2003). It also stands out in possessing the shortest 
period yet known among objects with black hole pri- 
maries (4.1 hr) and has exhibited only low luminosity 
hard state outbu rsts, r eaching just 1.2 x 10 36 ergs _1 
(IMcClintock et al.ll200l . 

The increasingly dominant paradigm for understand- 
ing the emission from XTE Jl 118+480, and from hard 
state BHXRTs in general, is that it (usually) involves 
an evaporated hot inner disk launchin g a compact jet 
l|Markoff et all 120011: lYuan et al.ll2005[h The jet is re- 
sponsible for not only radio emission, but also much 
of the IR, optical, and possibly even some UV cmis- 
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sion. The properties of the unusual UV, optical, 
and IR (UVOIR) variability also point to synchrotron 
emission rather than r eprocessing of X-r ays by a disk 
ijKanbach et~aT] I200H iHvnes et all Eoo.'i). In fact it 
has even been proposed that the X-ray emission seen 
in th e hard state could arise from synchrotron emis- 
sion ijMarkoff et al.l 120011) . although this is not widely 
accepted and may be ove rwhelmed by a la rger Comp- 
tonized X-ray component l)Yuan et alJl2005Tl . 

Part of the difficulty in bridging from the UVOIR re- 
gion to X-rays is that the former is likely a mixture 
of jet and disk emission, with the disk emission mask- 
ing the break from flat-spectrum to optically thin syn- 
chrotron. Consequently the position of the break, and 
slope of the optically thin component are not directly ob- 
servable, increasing the uncertainty in extrapolating to 
X-rays. The variability may provide the key to disentan- 
gling these components. To test this possibility, we have 
assembled serendipitous multiwavelength observations of 
XTE J1118+480 during the 2005 January outburst. The 
2005 outburst was a much shorter and somewhat weaker 
ev ent than that seen in 2000. It was discovered optically 
bv iZurita et alJ |2005a1. with the first high points seen 
on 2005 January 9. The ou tburst was also detected at X- 
ray and radio wavelength s ijRemillard et al.l2005t iPoolevI 
120051 iRuoen et alll2005lh The outburst f aded rapidly, 
reachin g near quiescence by late February ijZurita et alJ 
l2005b[) . For a discussio n of the outburst pr operties and 
longterm lightcurves see IZurita et alJ l)2006f) . 



2. OBSERVATIONS 

2.1. SARA 0.9 m Optical Observations 

XTEJ1118+480 was observed using the SARA 0.9- 
m telescope located at Kitt Peak National Observatory. 
Observations used an Apogee AP7p CCD camera with an 
R filter. Exposure times were 10 s, with approximately 
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7 s intervening dead-time. The raw data frames were 
bias-, dark-, and fiat field-corrected using standard IRAF 
routines. Once the data frames were calibrated, we ex- 
tracted time series aperture photometry using the exter- 
nal iraf package CCD_HSP written by Antonio Kanaan 
(U. Federal Santa Catarina, Brazil). CCD_HSP auto- 
mates the field alignment and photometry extraction for 
time-series CCD data. 

2.2. McDonald Observatory 2.1m Optical Observations 

XTE J1118 +480 was also observed usin g the Argos fast 
CCD camera ijNather fc Mukadamll2004|) on the McDon- 
ald Observatory 2.1m telescope. About 1.5 hrs of data 
were obtained in unfiltered white light. The data were 
taken as a continuous sequence of 1 s images with negligi- 
ble intervening dead-time. Conditions were mostly non- 
photometric with 1-2 arcsec seeing, so differential pho- 
tometry was performed relative to the same comparison 
star used in the IR. Data reduction employed a combi- 
nation of iraf routines to generate calibration files and 
then a custom IDL pipeline to apply calibrations and ex- 
tract photometry. Bias structure and dark current were 
subtracted using many dark exposures of the same dura- 
tion as the object frames. Residual time-dependent bias 
variations were removed using two partial bad columns 
which are not light sensitive. Unfortunately no flat fields 
were obtained in the unfiltered mode, so no sensitivity 
corrections were applied. Photometry was extracted us- 
ing standard aperture photometry techniques. 

2.3. McDonald Observatory 2.7 m Optical Observations 

Additional optical observations of XTE J1118+480 
were performed using the White Guider CCD camera on 
the McDonald Observatory 2.7m telescope. The data 
were taken as a continuous sequence of 3 s R band im- 
ages with about 2 s of intervening dead-time. Conditions 
appeared near-photometric with 1-2 arcsec seeing, so ab- 
solute photometry was performed to maximize the signal- 
to-noise ratio. Conditions deteriorated towards the end 
of the January 19 run, so these data were discarded. 
This approach was necessary as the comparison star used 
with Argos was not available with this instrument, and 
the brightest other comparison in the field was nearly 
six magnitudes fainter than XTE Jl 118+480 in outburst. 
Initial data reduction used standard iraf procedures to 
remove bias and flat-field the data. Aperture photome- 
try of XTE J1118+480 used a 5" aperture to minimize 
aperture losses. 

2.4. Kitt Peak National Observatory 2.1 m Infrared 
Observations 

Infrared photometry of XTE J1118+480 was obtained 
using the Simultaneous Quad Infrared Imaging Device 
(SQIID) on the 2.1m telescope at the Kitt Peak Na- 
tional Observatory. Details are given in Table^fe^ The 
camera records J, H , K Sl and L images simultaneously, 
although we found that the L data were of insufficient 
quality to be useful. JHK S data were all of good qual- 
ity, however, and obtained simultaneously with 1 s or 2 s 
exposures separated by ~ 54 s deadtime. The image was 
nodded back and forth on alternating exposures to fa- 
cilitate sky subtraction. The seeing was typically better 
than 1.3 arcsec in J. Data reduction employed a combi- 



nation of iraf 8 routines to generate calibration files and 
then a custom IDL pipeline to apply calibrations and 
extract photometry. 

Where possible, sky subtraction was performed by sub- 
tracting the average of immediately preceding and imme- 
diately following images. Sensitivity variations were cor- 
rected using an average of many sky flat images taken at 
twilight. We verified that these images acceptably flat- 
tened the sky background of target images before sky 
subtraction. Count rates were low enough that non- 
linearity was below 1 %, so no correction was applied. 

Because of the source brightness during outburst, only 
one usable comparison star was present in the field, and 
this was fainter than our target. We therefore used the 
comparison star only as a check of the photometric stabil- 
ity and did not perform differential photometry. Fortu- 
nately, conditions were mostly photometric, at least suffi- 
ciently so that transparency variations are much smaller 
than the intrinsic variability of the target. The major 
exception occured at the end of the third night, so these 
points were not used. 

Our absolute calibration is based on 2MASS photom- 
etry of the comparison star, 

2MASS J11180724+4803527 llCutri et alJ|200l . This 
star has J = 13.449 ± 0.024, H = 12.825 ± 0.028, and 
K s = 12.610 ± 0.022, where the uncertainties are domi- 
nated by statistical effects and so should be uncorrelated. 
We checked this cal ibration against A RNICA standards 
AS19-0 andAS19-2 ijHunt et alll 998) aft er transformin g 
the latter in the 2MASS JHK S system l|Caroenterl200Hh 
The standard yielded a consistent calibration, with com- 
parable uncertainties, so we retained the 2MASS calibra- 
tion for this work. Our average calibrated magnitudes for 
the target for each night are summarized in TableEl The 
ca librated magnitude s were converted to fluxes according 
to ICutri et alJ l)2003l) . After folding in uncertainties in 
conversion from magnitudes to fluxes, we estimate that 
the systematic uncertainties in our calibrated fluxes are 
then 2.8% in J, 3.2% in H and 2.8% in K s . 

2.5. Rossi X-ray Timing Explorer X-ray Observations 

XTE J1118+480 was intensively monitored in X-rays 
with the Rossi X-ray Timing Explorer (RXTE) from 
2005 January 13 to February 26. There are two 
narrow-field instruments aboard RXTE, the Propor- 
tional Counter Array (PCA) and the High-Energy X-ray 
Timing Experiment (HEXTE). Due to the faintness of 
the source throughout the outburst, only the more sen- 
sitive PCA produced useful data. For this work, we only 
used data taken in the Standardlb mode, which provides 
a time resolution of 1/8 sec and covers a nominal energy 
range of 2-60 keV. The PCA data were reduced with 
FTOOLS (version 5.2) that was distributed as a part of 
the software suite HEASOFT. 9 Briefly, for each obser- 
vation, we simulated background events with the appro- 
priate background model. The data were then filtered 
in the usual manner. 10 , which resulted in a list of Good 

8 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

9 see http://heasarc.gsfc.nasa.gov/docs/software/lheasoft 

10 see the online RXTE Cook book at 
http: / /heasarc. gsfc.nasa.gov/docs/xte/recipes/cook_book. html 
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TABLE 1 
Log of optical/IR observations 

Telescope & instrument UT date & time Exposure (s) 



SARA 0.9 m, AP7p 


2005 Jan 13 06:03:04- 


-13:50:50 


10 


SARA 0.9 m, AP7p 


2005 Jan 15 06:02:57- 


-13:53:07 


10 


McDonald 2.1m, Argos 


2005 Jan 13 07:42:07- 


-09:04:39 


1 


McDonald 2.7 m, White Guider 2005 Jan 16 07:51:06- 


-12:28:23 


3 


McDonald 2.7m, White Guider 2005 Jan 19 09:17:12- 


-10:56:48 


3 


KPNO 2.1m, SQIID 


2005 Jan 15 09:59:49- 


-14:00:29 


2 


KPNO 2.1m, SQIID 


2005 Jan 16 08:01:48- 


-13:52:09 


2 


KPNO 2.1m, SQIID 


2005 Jan 17 10:28:01- 


-13:51:59 


1 


KPNO 2.1m, SQIID 


2005 Jan 18 07:54:30- 


-13:54:12 


2 


KPNO 2.1m, SQIID 


2005 Jan 19 09:39:40- 


-13:53:08 


2 



TABLE 2 
Details of IR photometry 

UT date Number of images Average magnitude of XTE J1118+480 



All 


Good 


J 


H 


K B 


2005 Jan 15 221 


205 


12.92 


12.49 


11.97 


2005 Jan 16 360 


246 


12.96 


12.52 


12.05 


2005 Jan 17 210 


147 


13.03 


12.64 


12.10 


2005 Jan 18 380 


304 


13.06 


12.64 


12.14 


2005 Jan 19 270 


257 


13.10 


12.70 


12.19 



Time Intervals (GTIs). Using the GTIs, we proceeded to 
extract a light curve from the data and background files, 
respectively. The background-subtracted light curve was 
rcbinned to 1/4 sec for cross-correlation analysis. The 
only X-ray data used in this work are those with an over- 
lap with optical data. 

3. ORBITAL LIGHTCURVES 

The best data for searching for an orbital or super- 
hump modulation is the SARA photometry as this has 
the longest periods of coverage, and longer exposures 
suffer less scatter due to the large amplitude flickering 
present. A modulation on a period of about 4 hours is 
seen on both nights at about the same full-amplitude 
of ~ 2.0 percent (Fig. QJ. A period search including 
both nights of data finds several closely spaced aliases. 
One of them is at 0.17036 (24) days or 1.00 25(14) P nrb 
(assuming the period of iZurita et alJ 12*002'). This is 
consistent with previ ously reported super hump periods, 
e.g. 0.170529(6) days l|Uemura et al.ll2002fl but does not 
securely rule out an orbital modulation. Aliases at 
1.09 P rb and 1.20 P rb are statistically slightly preferred 
to that at the superhump period. The similarity to pre- 
vious observations, points to the solution closest to the 
orbital period being the correc t one, howe v er, as this 
is consistent with the finding of Ch ou et all l)2005[) of a 
very weak (0.02 mag) modulation close to the orbital (or 
superhump) period aft er January 18. I n con trast, the su- 
perhump identified bv lUemura et alJ l)2002ft throughout 
the 2000 outburst with exhibited a persistent modula- 
tion with an average full-a mplitude of about 7 %, so our 
observations and those of iChou et alJ l|2005D indicate a 
weaker superhump in this outburst. 

Where we have observations of at least one binary or- 
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Fig. 1. — Modulation in SARA photometry folded on preferred 
1.0025 P or b period, likely due to a superhump. 



bit in other data, we examined the lightcurves for evi- 
dence of a similar superhump modulation. These other 
lightcurves are less well suited as rapid variability is not 
averaged out. We find no evidence of any such mod- 
ulations in either the optical or IR data. The optical 
lightcurve from January 16 places a 2 a (95 %) upper 
limit on the full amplitude of a sinusoidal modulation at 
the orbital period (for any phasing) of 3.5 % of the flux. 
We attempted to analyze the IR lightcurves in the same 
way and, while no consistent modulation was apparent, 
the formal limits derived were much weaker. This is a 
consequence of the larger intrinsic variance due to flar- 
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Fig. 2. — Power density spectrum derived from fast photometry 
on 2005 January 13, in unfiltered white light. The dashed line is 
a single-Lorentzian component, the dotted line is the fitted white- 
noise level, and the solid line is the combined fit. Fits only used 
data above 0.01 Hz, where source variability is dominant. 
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Fig. 3. — Cross-correlation function between X-ray and 3 s time- 
resolution optical data from January 16. The inset shows an ex- 
panded view of the peak. 



ing, and the much lower duty cycle reducing the ability 
to average over many flares. Upper limits obtained were 
7-16% in J, 12-19% in H, and 17-24% in K s . 

4. POWER DENSITY SPECTRA 

Our best time-sampling is in the fast optical photom- 
etry from 2005 January 13, 1 s time-resolution with neg- 
ligible deadtime. In Fig. [21 we show the power-spectral 
density (PSD) of the differential lightcurve between the 
target and comparison star. Examining the PSDs of 
each star individually, the target exhibits excess power 
above about 0.01 Hz, even though it is brighter, while 
below that both stars have the same power, likely due 
to transparency variations and/or aperture losses. We 
are therefore confident that the power in the differen- 
tial PSD above 0.01 Hz should be due to real variations 
in the target, apart from white-noise which dominates 
at the very highest frequencies. We have shown the 
PSD in the now common vP v form and show a fit com- 
prising a single Lor entzian component plus white-noise 
l)Belloni et alJ f2002). The Lorentzian frequency param- 
eters are a characteristic frequency v maK = 0.13 Hz and 
central frequency v$ = 0.065 Hz. u max corresponds to 

the peak in vP v and is given by f max = vK) + ^ 2 wnere 
A is the half-width at half maximum of the Lorentzian. 
The difference from a zero-centered Lorentzian is small, 
corresponding to a low coherence Q = 0.11. The PSD 
derived is thus very similar to those seen during the 2000 
outburst, although the break frequency, z/ max , is higher 
than seen then, wh en it evolved from about 0.03-0.08 Hz 
(|Hvnes et alJl200^ . 

If the PSD did not change substantially between 2005 
January 13-19, and is similar in the IR and optical, then 
this gives us information about how well the IR obser- 
vations resolve variability. The single fitted Lorentzian 
component has 90% of its power on timescales 0.9-70 s. 
Thus most of the variability (> 90%) is on timescales 
shorter than the ~ 50 s sampling time, and hence con- 
secutive points are effectively independent. On the other 
hand, most of the variability is on timescales longer than 



the IR exposure time, hence the IR observations should 
not smooth out much of the variability and should sam- 
ple most of the dynamic range expected. 

5. X-RAY OPTICAL CORRELATIONS 

A small amount of overlapping X-ray and optical/IR 
data were obtained. The poor duty cycle of the IR data 
resulted in only a handful of points and no measurable 
correlation. A correlation was seen between X-ray and 
optical variability, based on a total of approximately 
llmin of overlap with McDonald 2.7m White Guider 
photometry. Fig. |21 shows the cross-correlation function 
(CCF) from 2005 January 16. Since the optical exposure 
times were quite long, this was calculated by averaging 
the X-ray data over the duration of each optical expo- 
sure (after applying the lag). The CCF is not of high 
quality, and the 3 s time-resolution (set by the exposure 
duration, not the cycle time) limits the precision of the 
information obtained, but a lagged correlation is clearly 
present. The centroid and width of the CCF pea k are 
comparable to those seen bv iKanbach et all l)2001|) dur- 
ing the 2000 outburst although the lower time-resolution 
of our data do not permit a detailed comparison. 

6. THE SPECTRAL ENERGY DISTRIBUTION (SED) OF 
THE IR NOISE 

6.1. Observational results 

Since J, H, and K s are observed simultaneously, each 
image-set provides an instantaneous IR SED. Compari- 
son of them allows us to study the SED of the IR vari- 
ability, even though we do not resolve individual flares in 
time. As a relatively large number of such image-sets was 
obtained, we applied quite cautious standards in filtering 
the data. We excluded all sets of images in which i) the 
seeing was greater than 1.7 arcsec, ii) the sky subtraction 
in any band was poor iii) the J sky was high in twilight, 
iv) conditions did not appear photometric (as measured 
by the comparison star), or v) the target coincided with 
either a hot pixel or an a-particle hit on the detector. Ta- 
bic lists the number of good image-sets remaining after 
applying these cuts to the sample. In Fig. 0] we show 
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examples of J — H and H — K s flux relations for the first 
and last nights and all nights with 2 s data combined. 
The cuts applied do appear to have cleanly removed bad 
data, leaving a well-defined correlations. We also show 
equivalent points for the comparison star. For the latter 
each night's fluxes have been divided by the median, so 
this is a measure of the scatter in fluxes within a night. 
If systematic effects, for example variable extinction due 
to cirrus, were affecting our data then we would expect 
the comparison star to exhibit correlated errors in dif- 
ferent bands. No such correlation is seen, suggesting the 
uncertainties are dominated by random statistical noise. 
The uncertainties in the comparison (which should be 
larger than those in the brighter target) are clearly much 
smaller than the observed variations in the target fluxes. 
The scatter about the correlation line in the target may 
be statistical, however. 

We performed the same analyses on the data from the 
third night, for which 1 s exposures were taken, to test 
if shorter exposures sampled a larger dynamic range of 
the variability. This does not obviously appear to be the 
case, which is consistent with most of the power being at 
lower frequencies (Section 0J . Since these data are not- 
icably noisier than the other nights, we have not shown 
them in the plots. The overall distribution of fluxes is 
intermediate between the preceding and following nights, 
with no obvious differences other than the higher noise 
level. 

Fig. 01 shows clear and repeatable correlations between 
bandpasses. Not only do the fluxes within a given 
night trace out a well defined relationship, but succes- 
sive nights appear to follow approximately the same cor- 
relation, with the overall brightness declining along the 
same line as the variations within a night. This points 
to a very consistent IR behavior and SED, and probably 
a single spectral component dominant. The consistency 
coupled with the linearity over a factor of two in flux 
suggest that the shape of the IR SED is not changing, 
only the overall flux level. 

We can use the correlations to estimate the shape of 
the variable component of the IR SED, since the slope of 
the correlation, for example dFx /dFjj, is a measure of 
the slope of the SED variations between H and K. This 
measure is independent of any non- varying zero-point, 
making this potentially a more sensitive discriminant be- 
tween models than the overall SED whic h may be a sum 
^of ^several^ spectra l components (e.g. iMarkoff et al]l200ll 

6.2. Power law models 

The simplest model for the IR SED is a power-law. 
We also tested broken power-laws (i.e. allowing a dif- 
ferent slope for J — H and H — K s ), but did not find a 
convincing pattern (see Table [3J). There is a tendency for 
the J — H slope to become flatter, and the H — K s one 
to begin steeper (corresponding to a concave spectrum) , 
but suspiciously the J — K s slope stays approximately 
constant, and tracks a single power-law fit to the full 
JHK S datasets. This suggests either that the H data 
are less reliable, or that adjacent bands simply do not 
provide adequate leverage in wavelength. We therefore 
see no compelling reason to believe the spectral slope 
really changes across the near-IR and proceed to fit a 
single power-law to the combined JHK S data for each 



Date 


ocjh 




ajK 


OtJHK 


x 2 / dofa 


Jan 15 


-0.75 ±0.12 


— 0.73 ± 0.11 


-0.78 ±0.07 


-0.77 ±0.06 


1.890 


Jan 16 


-0.90 ±0.08 


— 0.78 ± 0.09 


-0.85 ±0.04 


-0.83 ±0.04 


1.822 


Jan 17 


-0.30 ±0.16 


-1.04 ±0.15 


-0.67 ±0.08 


-0.70 ±0.07 


1.930 


Jan 18 


-0.50 ±0.09 


-0.98 ±0.08 


-0.74 ±0.04 


-0.76 ±0.04 


1.996 


Jan 19 


-0.48 ±0.07 


-1.02 ±0.08 


-0.74 ±0.06 


-0.74 ±0.05 


1.789 



a X 2 is for the joint JHK fit. 

TABLE 3 
Power law fit parameters 



night; the following discussion relates to this joint fit. 
Note that given three points in wavelength, such a bro- 
ken power-law model should be sufficient to completely 
describe the wavelength dependence, and this test can 
also be thought of as testing whether there is a signifi- 
cant and repeatable curvature of the spectrum. 

Free parameters in the unbroken power-law model are 
the assumed zero-point fluxes in each band (we adopt the 
mean), the normalization of the power-law for each im- 
age set, and the power-law slope. For n JHK image sets 
there are thus n + 4 free parameters and 2n — 4 degrees of 
freedom. The resulting fits typically yielded a \ 2 value of 
about 2 per degree of freedom, assuming only statistical 
errors. It is likely that the errors are underestimated in 
this way, as we did not perform differential photometry 
and hence aperture losses and small transparency vari- 
ations could contribute. To approximately correct for 
this, we adopt the expedient of rescaling the uncertain- 
ties (by a factor about y/2) to yield a minimum \ 2 °f 
1 per degree of freedom. We then estimate the single- 
parameter uncertainty in the slope, a (defined such that 
F v oc v a ), using the region encompassed by A\ 2 = 1. 
Our derived slopes are given in Table [21 and we also de- 
rive a weighted mean of the slopes of a = —0.78. The 
four nights are consistent with their weighted mean to 
within the quoted errors, and there is no overall trend to 
the values, suggesting a constant spectral index over the 
period observed. 

Given that the statistical uncertainty is small, and the 
night-to-night consistency is good, the dominant error 
will arise from the systematic uncertainty in absolute cal- 
ibration. Using our estimate of the absolute calibration 
uncertainties from Section l2~4l we estimate using a Monte 
Carlo calculation that the random uncertainty in the de- 
rived power-law index will be about ±0.07. This assumes 
errors in J, H, and K s are independent. If they are pos- 
itively correlated then the uncertainty in the power-law 
index will be reduced. Only if there is an anti-correlation 
between J and K s errors will we have underestimated 
the uncertainty in the power-law index, and this is un- 
likely. This calibration uncertainty dominates over the 
other terms so far discussed. The final problem that in- 
troduces a systematic bias rather than a random scatter 
is that the variability is redder than a stellar SED. We 
have treated the photometry as yielding monochromatic 
fluxes at the bandpass center, but the measurements ac- 
tually represent weighted averages over the bandpass. In 
pr actice t his eff ect is not large, however, so we neglect 
it. IGlassI (1999) estimate that the corrections to near-IR 
fluxes to convert them to the monochromatic flux at the 
bandpass center are only a few percent for power-laws of 
spectral index between —2 and ±2. We estimate with 



Hynes et al. 




2 4 6 a 10 12 140 2 4 6 B 10 12 140 2 4 6 S 10 12 14 
F (mJy) F (mJy) F (mJy) 



Fig. 4. — J vs. H and vs. K B fluxes for the first night, the last night, and all nights with 2 s data combined. Two clusters of points 
are present on each plot. The diagonally dispersed ones are XTE J1118+480, the circular clump are the comparison star. The dashed lines 
all show the same linear fit to the combined dataset (for that color), to provide a point of reference in comparing the first and last night. 
The cross indicates the run averaged fluxes as a point of reference. The joint fit shown here is also only intended as a reference point, and 
all analysis used night-by-night fits. 



some simplistic simulations that for our JHK S photom- 
etry, with the spectral slope derived, the resulting error 
in the spectral slope is likely to be less than 0.01, and 
much less than other uncertainties. We have therefore 
not attempted a more rigorous treatment. 

In summary, all of our observations are consistent with 
a single power-law of constant slope a = —0.78 ± 0.07, 
where the dominant uncertainty is in the absolute cali- 
bration of the photometry. 

6.3. Black body models 

A common paradigm for optical variability in X-ray 
binaries is that it arises in reprocessed X-rays. X-ray ir- 
radiation deposits energy at a modest optical depth in 
the atmosphere of the disk or companion star. This is 
then thermalized and emerges as optical and ultravio- 
let flux which would, simplistically, be expected to have 
something close to a black-body SED. This interpreta- 
tion was rejected based on detaile d analysis of the vari- 
ability during the 2000 outburst ijKanbach et alJ l200ll 
iHvnes et alJ |2003). Our results further support this, as 
the variability we see is clearly too red to arise in black 
body emission from the binary. Black body fits analagous 
to those discussed the preceding section yield tempera- 
tures ~ 1500 K, too low for this to be a plausible inter- 
pretation. 

6.4. Optically thin thermal models 



Although an optically thick thermal model is not plau- 
sible given the redness of the variability, optically thin 
thermal emission could still explain the spectral shape 
as the free-free emissivity increases with wavelength, re- 
s ulting in a redder SEP than that of a black body. 

iPearson et alJ l)2005tl derived analytic expressions for 
the time-dependent continuum spectra exhibited dur- 
ing flickering and flaring events. The additional flux 
was modeled as arising from a region of gas, with uni- 
form temperature and Gaussian density profile, expand- 
ing with a radial velocity proportional to the distance 
from the centre. Both free-free and bound-free emission 
mechanisms were considered and the results of fitting to 
observations consistently showed that the expansion ex- 
hibited isothermal evolution. In the optically thin limit, 
the observed flux is given by 



j? na2B ft\ 

F = ^ T0 (1) 

_ 7ra 2 B ( k\M 2 1 - erir \ 

= a„f (4) 
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Date 


Temperature 


X> 








Jan 15 


11 500 


1.888 


Jan 16 


10 600 


1.819 


Jan 17 


12 800 


1.940 


Jan 18 


11 700 


2.008 


Jan 19 


12 000 


1.797 




TABLE 4 





Derived temperature from the flickering occuring on 
each night for an optically thin thermal model. 



where a is the lengthscale of the Gaussian (= v2s.d.), B 
is the Planck function, tq is the optical depth through the 
centre of the expansion, d is the distance to the object, 
T is the temperature of the region and M its mass. The 
term k\ is a constant that depends on the composition 
and the emission mechanism (i.e. whether free-free or 
bound- free). a v and / encapsulate a number of terms 
from 1(3} , while separating the wavelength dependence 
into a v alone. 

Since the quantities contributing to / are fixed for a 
given observation of a particular flicker, the ratio of flux 
in each waveband is given by the ratio of the values for 
the ol v . Experience suggests that not only is T constant 
for a particular event but is also fairly consistent between 
flickers. The linear relationship between fluxes in the dif- 
ferent wavebands can thus be understood as a reflection 
of the linear relationship between the flux and the pa- 
rameter /. 

For each data triple, consisting of the flux in the JHK S 
wavebands, we can derive a best estimate instantaneous 
value for /. We can then fit a straight line to F(f) in 
each waveband simultaneously and extract a best fitting 
value for aj (or equivalently an or olk) which in turn is 
a simple function of T. In practice, rather than a single 
parameter, we also have to allow for a zero point offset 
in the flux in each waveband to account for constant 
or slowly varying contributions from other parts of the 
sytem. The derived values of T are given in table Elan d 
have x 2 /^ similar to the power-law fits in Section l6~2*l 

As a consistency check, we can insert the condition 
that material be optically thin (to < 1) into (JJJ. This 
gives us, 



Othin > 



2d 2 F 



ttB 
> 3 x 10 9 



(5) 
(6) 



using d = 1.8 kpc, T ~ 12 000 K and F K ~ 20 mJy. 
Unfortunately then, while the thermal models do manage 
to reproduce the data equally as well as a power law, 
they require an emiting region comparable in size to the 
entire binary (a ~ 2 x 10 9 m). It is hard to envisage 
how this could arise from material within the primary's 
Roche Lobe. 

This difficulty suggests that, despite the attractiveness 
of the fit, the emission mechanism, in this case, is not 
thermal in origin, and hence that our original power-law 
fits are probably most meaningful. 

7. DISCUSSION 

The properties of the variability we observe are very 
similar to those seen during the 2000 outburst, for ex- 



ample the prompt correlation with X-rays, the shape of 
the PSD, and t he increasing ya riability at longer wave- 
length s. Both iKanbach et al.l l)2001|) and iHvnes et alJ 
( 2003) concluded that synchrotron emission was the most 
likely origin of the variability. Our multicolor obser- 
vations provide strong support for this interpretation. 
We have argued that a thermal model, whether optically 
thick or thin, cannot adequately explain the combination 
of spectral shape, luminosity, and variability timescales 
observed. The inferred IR spectral slope, a ~ —0.78 is, 
on the other hand, quite appropriate for optically thin 
synchrotron emission. 

In the context of the jet model for the broad-band SEP 
of XTEJ1118+480 presented bv iMarkoff et alJ ^2001^ . 
this implies that the IR variability, and probably also 
the optical and some UV, originate from the "optically 
thin post-shock jet" component . This may i ndicat e a 
difference from the models of IMarkoff et alJ <)2001[) in 
which the near-IR exhibits a flat spectrum, as this implies 
that the break between optically thin and partially self- 
absorbed synchrotron occurs at longer wavelengths than 
the near-IR. The slope we derive, however, a ~ —0.78, 
appears in good agreement with these models. Our ob- 
servations therefore provide new tests of these models, 
allowing isolation of the optically thin jet emission from 
c ontamination by th e disk. 

lYuan et aTl l)2005l) consider a coupled accretion-jet 
model. In this model, the IR emission is mostly due 
to radiation from the jets, although the contribution 
from hot accretion flows might not be negligible. The 
jet contribution is associated with synchrotron emission 
from electrons that are accelerated by internal shocks 
in the jets. The spectral distribution of the electrons 
is assumed to be of power-law shape, N{E) oc E~ p for 
£min < E < E max , where p = 2.24. For optically thin 
synchrotron emission one then expects that the spectrum 
of the radiation is also of power-law shape, for the most 
part, with a spectral slope a = (1 — p)/2 = —0.62, which 
is already rather close to that which we observe. More- 
over, for XTE J1118+480, the IR emission may be asso- 
ciated with electrons near i? max , so its spectrum devia tes 
from the power law and is steeper ijYuan et al.ll2005]) . 

8. CONCLUSIONS 

We have performed optical, IR, and X-ray observa- 
tions of rapid variability in XTE Jl 118+480 during its 
2005 outburst. Many characteristics are similar to the 
2000 outburst, although there are a few key differences. 
Superhumps, if present, are weaker than in 2000. The 
variability also seems concentrated to higher frequencies. 
Our major novel result is simultaneous J, H, K s photom- 
etry allowing us to isolate the SED of the IR variability. 
We find this is red, and can be well fitted by a power- 
law, F u oc v a , where a = -0.78 ± 0.07. This result is 
consistent with optically thin synchrotron emission but 
hard to explain with thermal emission. We consider this 
to be strong evidence in favor of the interpretation of the 
variability as arising in synchrotron emission, most likely 
from a jet. Unlike attem pts to model the average spectral 
energy distribution (e.g. IMarkoff et alJl200H [Yuan et all 
2005), the variability isolates the synchrotron emission 
from the disk emission allowing us to measure of the 
slope of the optically thin component directly and show 
that the spectral break to self-absorbed emission must 
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occur at longer wavelengths than the near-IR, at least at 
the time our observations were made. 
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